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^Rapid  melting  and  solidification  of  a  semi-infinite  substrate 
subjected  to  a  high  intensity  heat  flux  over  a  circular  region  on 
its  bounding  surface  moving  with  a  constant  velocity  is  considered. 
General  expressions  are  developed  for  the  coefficients  in  the  fin¬ 
ite  difference  equation  governing  the  heat  transfer  in  moving 
orthogonal  cuvilinear  coordinate  systems.  These  expressions  are 
reduced  to  their  specific  forms  in  terms  of  dimensionless  nodal 
temperature  and  enthalpy  for  a  moving  oblate  spheroidal  coordinate 
system.  Quasi-steady  state  conditions  are  assumed  and  the  thermal 
properties  of  the  substrate  in  the  liquid  and  solid  phase  are  con¬ 
sidered  constant  and  equal.  It  is  also  assumed  that  the  substrate, 
pure  aluminum  used  as  example,  melts  and  solidifies  at  a  single 
temperature.  Temperature  distributions  in  the  molten  region  and 
the  adjacent  heat  affected  zone  are  computed  along  with  the  liquid- 
solid  interface  shape,  its  velocity  and  other  important  solidifica¬ 
tion  variables.  Both  uniform  and  Gaussian  heat  flux  distributions 


within  the  circular  region  are  considered.  The  results  are  present¬ 
ed  in  their  most  general  form  -  in  terms  of  dimensionless  numbers 
v< hen  possible.  Specific  criteria  for  the  melting  of  the  substrate 
are  established,  ,/jlt  is  shown  that  the  three  variables,  absorbed 
heat  flux  the  Radius  of  the  circular  region  a_  and  the  velocity 
of  the  moving  heat  flux  U,  could  be  combined  into  two  independent 
variables.  That  is,  the  dimensionless  temperature  distribution  in 
the  metal  pool  and  the  solid  substrate  remain  the  same  as  long  as 
the  products  <ja  and  U_a  or  U/q  are  kept  con-stant.  The  effect  of 
these  variables  on  cooling  rate  in  the  liquid  and  the  ratio  of  temp¬ 
erature  gradient  to  growth  rate  at  the  solid-liquid  interface  are 
discussed  using  the  aluminum  substrate  as  example. 
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I.  INTRODUCTION 


The  availability  of  high  power  directed  energy  sources  such 
as  the  electron  beam  and  different  types  of  lasers  has  led  to 
the  development  of  a  number  of  new  materials  processing  techniques 
which  exploit  the  unique  characteristics  of  these  sources.  One 
such  process,  the  rapid  surface  layer  melting  and  subsequent  sol¬ 
idification  of  metallic  and  semi-conductor  substrates,  appears 
to  have  tremendous  potential  applications.  In  two  recent  papers 
(1,2)  we  addressed  the  one  and  two  dimensional  transient  heat 
flow  problems  during  rapid  melting  and  solidification  of  the  sur¬ 
face  of  a  semi -i nf i ni te  substrate  subjected  to  a  high  intensity 
stationary  heat  flux  on  its  bounding  surface.  In  the  present 
investigation  we  extend  the  earlier  findings  to  three-dimensional 
heat  flow  on  the  surface  of  a  semi-infinite  solid  subjected  to  a 
moving  heat  flux.  It  is  anticipated  that  the  equations  and  sol¬ 
ution  method  developed  would  be  equally  applicable  to  other  met¬ 
allurgical  processes  such  as  arc  welding. 

In  general,  most  experiments  with  a  directed  energy  source, 
such  as  the  continuous  wave  CO2  laser,  involve  scanning  of  the 
source  over  the  surface  of  the  substrate.  Analytical  solutions 
to  simple  moving  heat  source  problems  have  previously  been  con¬ 
sidered  by  Rosenthal  (3).  His  analysis  is  for  a  solid  substrate 
which  does  not  undergo  a  phase  change.  It  is  based  on  the  notion 
that  if  the  dimensions  of  the  substrate  are  large  with  respect 
to  the  moving  source,  then  the  system  approaches  a  quasi-steady 
state;  steady  state  prevails  from  the  standpoint  of  an  observe** 
located  in  and  travelling  with  the  source.  The  analytical  solutions 
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of  Rosenthal  (3)  have  been  extensively  used  in  metallurgical 
processes  such  as  welding  and  surface  hardening.  However,  these 
solutions  are  only  accurate  at  large  distances  from  the  source 
and  can  not  address  the  complex  problem  of  melting  and  solidifi¬ 
cation  which  Is  the  subject  of  this  investigation. 

In  this  paper  we  extend  the  mathematical  technique  developed 
and  used  in  the  previous  two-dimensional  transient  heat  flow 
problem  (2).  The  oblate  spheroidal  coordinate  system  is  used 
again,  however,  the  mathematical  expressions  and  computer  method¬ 
ology  developed  assume  the  existence  of  quasi-steady  state  while 
the  coordinate  system  is  in  motion. 
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II.  PROBLEM  STATEMENT  AND  SOLUTION  APPROACH 

We  consider  a  high  intensity  heat  flux  over  a  circular 
region  on  the  bounding  surface  of  a  semi-infinite  solid  moving 
with  a  constant  velocity,  U,  in  the  y-direction  in  cartesian 
coordinate  system.  Figure  1.  The  absorbed  heat  flux  is  high 
enough  to  cause  melting  and  subsequent  solidification  of  the  sur¬ 
face  layer.  Temperature  profiles  in  the  molten  region  and  the 
adjacent  heat  affected  zone,  as  well  as  the  important  melting  and 
solidification  variables  of  the  surface  layer  are  to  be  determined. 
The  analysis  is  based  on  the  assumption  that  a  quasi-steady  state 
is  established,  which  is  to  say  that  the  system  appears  to  be 
steady  state  as  viewed  by  an  observer  located  at  the  center  of 
the  circular  region  and  travelling  with  the  heat  source.  We  thus 
transfer  the  coordinate  system  from  the  semi-infinite  solid  to  the 
center  of  the  heat  source.  The  surface  outside  the  heated  region 
is  considered  adiabatic.  The  thermal  properties  of  the  solid  and 
the  liquid  phases  are  considered  to  be  constant  and  equal  to  one 
another.  Finally,  it  is  assumed  that  the  workpiece  melts  and  sol¬ 
idifies  at  a  single  temperature. 

The  generalized  expressions  previously  derived  (4)  for 
the  determination  of  the  coefficients  in  the  finite  difference 
equations  governing  the  stationary  heat  transfer  problem  within 
discretized  spatial  domains  are  extended  to  account  for  the 
motion  of  orthogonal  curvilinear  coordinate  systems.  These  ex¬ 
pressions  are  then  reduced  to  their  specific  forms  for  a  moving 
oblate  spheroidal  coordinate  system  which  is  a  "more  natural" 
coordinate  system  for  this  problem  geometry.  The  finite 


4 


difference  equations  are  rewritten  In  terms  of  dimensionless 
nodal  enthalpy  and  temperature.  In  a  manner  similar  to  that 
previously  described  (2,6),  to  permit  numerical  solution  of 
the  multidimensional,  discrete  temperature -phase  change  prob¬ 
lem.  In  this  way,  both  dimensionless  temperature  and  enthalpy 
are  used  to  formulate  a  single  energy  conservation  equation  for 
each  discretized  spatial  domain  regardless  of  whether  it  is  in 
the  solid  state,  the  liquid  state  or  contains  the  liquid-solid 
Interface.  Finally,  the  quasi-steady  state  temperature  distri¬ 
butions  in  the  molten  region  and  the  adjacent  heat  affected 
zone  are  computed  along  with  the  liquid-solid  interface  shape, 
its  velocity  and  other  important  melting  and  solidification  var- 
iabl es . 
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V 

III.  MATHEMATICAL  DESCRIPTION 


The  generalized  form  of  the  heat  conduction  equation  in 
stationary  orthogonal  curvilinear  coordinate  system  (u^,  Ug»  u.j) 
has  previously  been  derived.  For  a  volume  element  moving  in 
space  with  velocities  v.| ,  V2  and  v.j  this  expression  can  be  ex¬ 
panded  to  include  conduction  of  heat  in  and  out  of  the  volume 
element  due  to  motion:* 


3U,  l  7?  3U,J  3U2  L  au2  J  3-3 


(1) 
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Where  the  scalar  factors  (metric  coefficients)  relating  the 
curvilinear  coordinate  system  to  the  cartesian  system  are  those 
previously  defined  (  5): 


The  velocities  in  the  two  coordinate  systems  are  related  by 


(  3) 


,  =  _L  (i2L_)v  +  —  (-$JL)v  +  -r—  (JLi.)v  • 
i  h.  '3ui;vx  h.  l3u.;vy  h.  ^3u.'vz* 

1  =  1,2,3 


(4) 


♦All  the  terms  in  the  equations  are  defined  in  the  Nomenclature. 
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The  arc  lengths,  areas  and  the  volume  of  the  element  in  the 
curvilinear  coordinates  are  related  to  the  cartesian  by: 


ds.  =  h.dui 


i  =  1  ,2,3 


dA.  =  h.h.  du.du.5  i.j»k  =  1*2,3  (6) 

1  3  K  J  K 

dV  8  h  dUj  du2  du3  (7) 

Finally,  P  in  equation  (1)  denotes  source  strength  per  unit 
volume. 

Taylor  series  expansion  of  equation  (1)  about  an  arbitrary 
discretized  domain  in  space  centered  at  a  node  at  (i,j’,k)  can 
now  be  carried  out  to  put  the  equation  in  finite  difference  form. 
The  subscripts  i,j  and  k  indicate  the  finite  discretization  of 
space  in  the  u^  ,  u2  and  u3  directions,  respectively.  The  change 
in  respective  coordinate  values  between  successive  nodes  are 

3 

AOj  Au2  and  Au^.  If  terms  of  the  order  of  Au  and  higher  are 
neglected,  the  following  expressions  are  obtained. 


r  -H- 1 
L  .  T  an  J 


Ti*i 


i.j.k 


■k  *  Ti-i 
(Aut) 


j ,k  '2Ti  ,k 


_3_  klh  rTH-1  ,3  ,k  '  Ti-l,j,k1 

lu,  “IT  *■  2Au,  J 

1  "l  i.j.k  1 


p.h  -  [P.h])>jik 


>  (pC_7).  i  k  ”  (P^n^i  i  |< 

Ti  <PcphT)  =  h  - £-UbjL 3 

1,j,k 
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where  the  superscript  refers  to  the  previous  time  level. 


-  h  3T  ,  ~  _h_  u 
pCp  h1  V1  3u1  pCp  h]  V1  * 


rTi^l,j,k  ~  Ti-l.J.k 
t  2Au.  J 


i.j.k  (11) 

Similar  expressions  to  (8)  -  (11)  are  then  written  for  u2 
and  Uj.  Substitution  of  all  the  finite  difference  approximations 
in  equation  (1)  yields  the  coefficients  to  the  various  nodal 

temperatures : 

r !  k-,  h  ,  „  k,  h  .  pCnhv,  i{ 

I  /  r  I 


f  Ki 

IhT 


h^  Au^ 


J _ 3_  /  kih  i  _J_ 

-Au .j  9u-|  1  ^T~  ‘  '  2Au 


1  J 


li  »j  »k 


h2Au2 


_1 _ 3_  /  ^2_h  x  1  /  pCphv2  x 

2Au2  3u2  (  }'  ^7  (  ) 


Ti  .j- 1  ,k 


i.j.k 


,  r  k3h  1  9  ,  k3h  x  1  ,pCphV3 

"  2Au3  3u3  h^~  2Au3  *  h3 


Pn  3.' 

^  }\ 


i  .j  ,k-l 


|i .  j.k 


,  r  Kin  ,  i  3  ,ki\  ,  i  fpCphvi 

^  ^  *=r  ‘if  (  ^ 


li  ,j  ,k 


r  k aH  i  ^  k~h  -j  pC  hV/ 

+  ^  V  +  ***  (_V 


i+1  ,j  ,k 


i .  j  +  1  ,k 


i.j.k 


+  [  k3h  1  _9_  .  ,  _J__  , 

7T~?  2auT  9u„  [  u2  J  2Au„  1 


pc  hv_ 

D  3 


h?U3 


3  3  h. 


Ti,j.k+1 


i.j.k 


k-|h  k2h  k3h 

2  2  +  ~2  2  +  2  2 
h*Au^  h2iu2  h^AUg 


hPC_ 

-  <TF*) 

i.j.k 

i.j.k 

Ti.j.k 


+  P.h  +  ( 


i.j.k 
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Application  of  the  general  expression  (12)  to  a  moving 
oblate  spheroidal  coordinate  system  located  at  the  center  of 
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the  heat  flux  and  travelling  with  the  source  is  now  considered. 
The  cartesian  and  the  oblate  spheroidal  coordinate  systems  are 
shown  in  Figure  l.  The  heat  flux  is  applied  over  a  circular 
region  of  radius  a^  in  the  x-y  plane  and  is  travelling  in  the 
y-direction  with  constant  velocity  U. 


(13) 


The  thermal  properties  are  considered  to  be  uniform  iso¬ 
tropic  and  equal  for  the  solid  and  liquid  phase.  Quasi-steady 
state  conditions  are  assumed  -  terms  involving  time  in  equation  (12) 
are  zero. 

The  relationships  between  the  cartesian  and  the  oblate 
spheroidal  coordinates  are: 

x  =  a  coshn  sin?  c o s 4> 
y  *  a  coshn  sin?  sin<f>  •  (14) 

z  =  a  sinhn  cos? 

The  control  volume  is  centered  about  an  arbitrary  point 
(i»  j»  k)  in  space  and  is  moving  with  a  velocity  U  in  the  pos¬ 
itive  y-direction.  i,  j  and  k  indicate  the  finite  discretiza¬ 
tion  of  space  in  the  n,£  and  $  directions,  respectively. 

u-j  -  nJ  u2  =  ?,  u3  - 


(15) 


\ 
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The  scalar  factors  (metric  coefficients)  for  the  oblate 
spheroidal  coordinates  determined  from  equations  (2)  and  (14)  are: 


h  =  hr 
n  € 


'cosh“n  -  sin' 


h^  =  a  coshn  sin?  (1 

h  =  h^.h^.h^  =  a(cosh2n  -  sin2?)  coshn  sin? 


The  interrelationship  between  vector  components  (velocitiss) 
in  the  two  coordinate  systems  from  equation  (4)  are: 

vn  =  (f£>  vx  *■  ^  vy  *  (§S>  vz 

^  (fc>  vx  *  h^T  vy  +  vz  07) 


i  (<**_)  v  +  _L  (<hL)  v  +  _L  ($-?.)  v 

h^  vx  +  h^  vy  h^ 


Substituting  equations  (14),  (15)  and  (16)  in  equation  (17) 
and  performing  the  indicated  operations  yields: 


=  U  cos0 

The  scalar  factors  (16)  and  the  velocity  components  (18)  are 
now  used  in  the  general  expressions  for  the  finite  difference 
coefficients  (12).  Once  the  appropriate  operations  are  carried  out 
expression  (12)  can  be  put  in  a  more  useful  form  for  the  problem 
at  hand. 
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Cs  Ti,j,k 


-  C,  T, 


1  '  i - 1 ,  j  » k  +  C2  1 i  + 1 , j , k  T  u3  1  i ,  j - 1 ,  k  +  u4  Ti,j  +  l,k  + 


+  ^ i  i  L  +  Ci.  T  .  . 


Cc  T 


+  Cf  T,  . 


5  ' i » j , k-1  T  ^6  1 i » j , k+1  "  C7  T  i  - 1 , j  ,  k  +  L8  'i  +  l,j,k  ~ 


+  C0  T, 


Co  T,  , 


"9  'i.j-l.k  +  C1 0  Ti ,  j  +  1 ,k  -  U11  'i,j,k-l  +  u12  Ti,j,k+1 
Pa2(cosh2q.  -  sin2?.)  .  . 


T .  , 


+  C„  T .  , 


The  coefficients  in  equation  (19)  are: 


C 1  , 2  =  An7  +  2Aq  tanhni 


cot? . 


T  m 

( cosh2^.  -  sin2?.) 
cosh2n.  sin2?. 

*  J 


alipC 


(20) 


a 

h  -  c8  *  HT&  s1"hr,i  5i"5j  s1n<t,k 

/ 

aU  pC 

c9  =  C10  '  21T4C  c°shni  cos5j  s1"*k 

aOpC  (cosh2ni  -  sin2?.) 

c  -  c  =  P  -  - - J  co  s  <j> 

11  12  2k  coshrii  sin?  .  k 

6  r  1  7  1  (coshzn,  -  sin2?.) 

Cs  -J,  Cn  ■  *2  ^  +  &  ♦  jir  - — 3  7  — J  3 

n-l  s  v  coshzn.  sin2?. 

^  J 


In  the  above  each  coefficient  has  been  divided  by  (ak  coshq. 
si n? j  AnA?A$). 

The  finite  difference  representation  of  the  heat  conduction 
equation  can  now  be  put  in  its  equivalent  enthalpy  form  using  a 
modified  version  of  previous  notation  (  2  )  for  dimensionless  nodal 


A  discretized  space  volume  containing  the  liquid-solid  inter¬ 
face  is  at  the  melting  point  of  the  material: 

0  «  ¥  $  1 .0  and  6=0  (24) 

The  value  of  is  equal  to  the  weight  fraction  of  the  element 
which  is  in  the  liquid  state,  f^. 
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The  appropriate  forms  of  equation  (21)  for  a  discretized 
space  domain  in  the  solid  or  liquid  region,  equations  (22)  or  (23) 
are  substituted  into  equation  (19).  After  some  manipulation  and 
multiplication  by  different  factors  for  the  following  common  equa¬ 
tion  is  obtained  for  both  the  solid  and  the  liquid. 


C  tb .  .  . 
s  M  ,j  ,k 


C1  ^i-U.k  +  C2  ^  i  + 1 ,  j ,  k  +  C3  ’h.j-l.k  +  C4  *i,j+l,k 
+  C5  6i,j,k-l  +  C6  6i,j,k+l  '  C7  Vl,j,k 
+  C8  9i+l,j,k  "  C9  8i,j-l,k  +  C1 0  9 i , j  + 1 , k 

(25) 

•  C11  9i  »3  «  k-1  +  C12  6i,j,k+l 
+  Pa2(cosh2ni  -  sin2?.)  C 

k  AHsr 


The  coefficients  to  and  C$  are  those  defined  in  ex¬ 
pressions  (20). 

An  alternate  approach  to  the  general  formulation  of  the  prob¬ 
lem  developed  here,  which  resulted  in  equations  (19)  and  (20) ,  is  to 
apply  an  energy  balance  to  an  arbitrary  moving  control  volume  of 
finite  size  with  a  total  source  strength  PAV  centered  about  node 
i,j,k  in  oblate  spheroidal  coordinates.  This  approach,  described 
in  detail  in  the  Appendix,  is  very  useful  because;  (a)  it  verifies 
the  general  formulation  of  equation  (1)  developed  in  this  investi¬ 
gation,  and  (b)  it  renders  improved  physical  interpretations  of  the 
various  coefficients  derived  for  equations  (19)  and  (20)  and  the 
boundary  conditions  described  in  the  next  section. 


IV.  BOUNDARY  CONDITIONS  AND  SOLUTION  OF 


THE  FINITE  DIFFERENCE  EQUATIONS 


The  boundary  conditions  are  derived  based  on  the  following 
assumpti ons . 

The  absorbed  heat  flux  in  the  circular  region  of  the  bound¬ 
ing  surface  (z  =  0,  x-y  plane)  is  in  general  a  function 
of  distance  and  time,  q=q(r,t)  where  r  =  /x *  +  y 7 .  In  this 
paper  the  problem  is  solved  for  the  two  special  cases  of 
uniform  and  Gaussian  distributions  of  the  absorbed  heat 
flux  within  the  circular  region.  Equation  (25)  is  then 
subject  to  the  following  boundary  conditions. 


I.  n  =  0  ,0  «  £  «  tt/2  ,  -  |  <  <J>  4  jrt  =  0  (26) 

The  first  node  at  this  boundary  is  at  i  s  1,  ^  8  4p-,  there 
fore,  the  surface  located  at  i  -  1/2  is  coincident  with  this  boun 
dary  of  symmetry: 


Q  i  - 1  /  2 ,  j  ,  k 


(27) 


Similarly,  the  top  surface  is  parallel  to  the  moving  direc¬ 
tion,  hence: 

C?  -  0  (28) 


In  general,  the  heat  flux  in  the  circular  region  on  the  bound¬ 
ing  surface  is  incorporated  in  the  term  involving  P  in  equation 
(25). 


P 


rate  of  heat  generated 
unit  volume 


€,+A?/2  $+A<|>/2 

a2/  J  q(5)  sin(2£)  d £./  d<f> 

7  E-AS/2  <j>-A<j>/2 

J 


(29) 


I 

14 

where  AV  is  the  volume  of  the  discretized  space  domain. 

AV  =  hAnA?A4>  =  a3(cosh2ni  -  sin2? Jcoshn.  sin?.AnA?A$  (30) 

I  J  I  J 

For  case  (a)  q  =  constant 

P  is  calculated  from  equations  (29)  and  (30).  Substitution  of  this 
finding  in  the  last  term  of  equation  (25)  leads  to  the  following 
expression : 


a  q  cos?,  sin (A?) 

J 

k  coshri^.AnA? 


(31) 


For  case  (b)  q  =  qQ  e"Z  sin 

and  the  last  term  of  equation  (25)  becomes: 


(32) 


coshn.  sin?.  AnA? 

1  J 


AH 


s£ 


(33) 


II.  n>0,0<?^Tr/2,  $  =  -  tr/2 ,  |I  *  o,  P  =  0 

(the  negative  y  portion  of  the  y-z  plane,  x  =  0)  (34) 

Because  of  problem  symmetry  about  the  y-z  plane  the  plane  defined 
by  <j>=-tt/2  will  represnet  a  zero  flux  boundary.  Since  the  first 
node  at  this  boundary  is  at  k  8  1,  ^  =  A4>/2,  then  the  surface 
located  at  k  -  1/2  is  coincident  with  this  boundary  of  symmetry: 


Q  i ,  j  ,  k  - 1  /  2 


k=l 


(35) 


Similarly,  this  surface  is  parallel  to  the  moving 
d i recti  on  ;  hence : 


C 


11 


0 


(36) 
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III.  n  >  o.  o  «  e  <  J.  ♦  ■  |I  -  0.  P  -  0  (37) 

(the  positive  y  portion  of  the  y-z  plane,  x=0) 

By  using  a  similar  reasoning  to  that  given  above  we  must  set 


Cfi  =  0 


C1 2  =  0 


n  -*■  °°,  0  <  £  <  w/2,  -tt/2  <  4>  <  tt/2  ,  P  =  0  (39) 

Far  away  from  the  circular  region  (for  the  problem  at  hand. 


n  10): 


6  i » j » k  "  ^  i  •  J  ,  k 


-CP  <TM  -  V 


V.  n  >  0,  5=0,  P  =  0  (41) 

(along  the  z-axis,  x  =  0,  y  =  0) 

The  area  of  the  surface  located  at  j  -  1/2  along  this 


boundary  is  zero 


Cg  =  0 


VI.  n  >  0,  £  =  tt/2,  -tt/2  <  <j>  <  tt/2,  P  =  0  (43) 

The  surface  on  the  x-y  plane  outside  the  circular  region 
is  adiabatic: 


Qi » j+l/2,k 


j=Max 


=  0  ,  c4  =  0 


The  top  surface  is  parallel  to  the  moving  direction,  hence: 

C10  -  0  (45) 


The  system  of  quasi-steady  state  algebraic  equations, 
equation  (25),  in  the  moving  oblate  spheroidal  coordinates 
were  solved  using  an  iterative  method.  The  computer  logic 


presented  below  closely  follows  that  previously  described 
for  the  two-dimensional  transient  heat  flow  problem  for 
stationary  heat  flux  applied  in  the  circular  region  on  the 
surface  of  a  semi-infinite  solid  (2).  However,  the  quasi¬ 
steady  state  nature  of  the  problem  has  eliminated  time  deriv¬ 
atives  from  the  heat  flow  equation. 

The  solution  is  started  by  initially  assigning  a  temper¬ 
ature  of  Tq  to  the  semi -i nf in i te  solid.  Then,  by  using  equa¬ 
tions  (22)  to  (24)  the  left  hand  side,  L.H.S.  of  equation  (25) 
is  calculated  using  the  boundary  conditions  by  repeated  point  iter¬ 
ation  throughout  the  mesh  in  a  definite  order  a  number  of  times 
until  the  convergence  criteria  is  met.  The  following  logic  is 
used  in  the  sequence  of  calculations. 

If  i  <  0  the  element  is  in  the  solid  and  right  hand 
side,  R.H.S.,  of  equation  (25)  is  less  than  zero. 

.  _  R.H.S. 

i . j  "  Cs  (46) 

As  the  calculation  is  repeated  for  the  next  nodal  point, 
the  value  of  0  in  the  previous  nodal  point  in  the  mesh  is  set 
equal  to  that  calculated  from  equation  (46).  On  the  other  hand, 
i  f  0  «  ip  .  .  .  $  1  . 0  the  element  contains  the  liquid  solid  inter- 
face  and  the  value  calculated  from  equation  (46)  gives  the 
fraction  of  liquid  in  the  volume  element.  The  value  of  e  for 
this  nodal  point  is  set  equal  to  zero  in  the  next  iteration  step 
Finally,  if  ij».  .  .  >  0  the  element  is  in  the  superheated  liq- 

1  9  J  9  K 

uid  region  and  the  values  of  <p  and  of  0  for  this  nodal  point 
become  that  given  by  equation  (23). 

The  convergence  criteria  is  tested  by  comparing  the  new 


value  of  with  the  old  guess  value: 
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V.  RESULTS  AMD  DISCUSSION 

The  equations  and  the  computer  logic  developed  were  used  to 
calculate  the  quasi-steady  state  temperature  distribution  in  an 
al uminum* substrate  subjected  to  both  uniform  and  Gaussian  moving 
heat  flux  distributions.  The  results  are  presented  in  their  most 
general  form,  in  terms  of  dimensionless  numbers  when  possible, 
in  order  to  establish  general  trends  between  the  process  vari¬ 
ables  and  the  important  melting  and  solidification  parameters. 

The  sequence  of  the  presentation  is  as  follows.  First,  the 
steady  state  temperature  distributions  due  to  a  stationary  heat 
flux  acting  over  a  circular  region  of  radius  a^  are  discussed  and 
the  results  are  compared  with  the  transient  heat  flow  calculations 
of  the  previous  paper  (2).  It  is  shown  that  the  criteria  developed 
earlier,  between  the  product  of  the  absorbed  heat  flux  and  the 
radius  of  the  circular  region  c£  and  the  steady  state  temperature 
at  the  center  of  the  circular  region,  are  equally  applicable  to 
the  problem  on  hand.  A  significant  departure  from  the  earlier  (2) 
calculations  is  the  assumption  that  the  conductivities  of  the 
liquid  and  the  solid  phases  are  equal.  The  effect  of  this  assump¬ 
tion  on  the  steady  state  temperature  distributions  is  discussed. 
Second,  the  effect  of  moving  the  heat  flux  in  the  y-direction  with 
a  dimensionless  velocity  Ua/2ot  on  the  temperature  distributions 
and  the  solidification  parameters  are  discussed  in  detail  for  both 
uniform  and  Gaussian  heat  flux  distributions. 

*The  properties  of  aluminum  used  in  the  calculations  are  listed 
in  Table  I . 


19 


1 .  Steady  State  Temperature  Distributions 
Stationary  Heat  Flux 

Figure  3  shows  a  general  plot  of  the  data  obtained  in  the  previous 
study  (2).  The  curve  associated  with  the  vertical  axis  on  the 
right  side  of  this  Figure  shows  that  there  is  a  minimum  product 
of  3^  required  if  the  center  of  the  circular  region  on  the  surface 
of  the  substrate  is  to  reach  a  given  temperature,  e.g.  the  vapor¬ 
ization  temperature  of  the  substrate.  That  is,  for  very  small 
values  of  a/2/at ,  long  interaction  times,  the  temperature  at  this 
location  approaches  its  maximum  steady  state  value.  Again,  the 
term  AHsjyCp  in  the  numerator  on  the  right  hand  vertical  axis  of 
Figure  3  denotes  the  equivalent  temperature  change  for  the  melt¬ 
ing  of  the  substrate. 

For  the  aluminum  substrate,  the  minimum  values  of  £a  c  1.45  x 
5  5 

10  W/m  and  =  2.3  x  10  W/m  are  deduced  from  Figure  3  for  a 

solid  surface  temperature  T(0,0,0)  =  TH  and  for  the  initiation  of 

surface  melting,  respectively.  These  values  are  identical  to 

those  of  the  previous  calculations  (2).  On  the  other  hand,  the 

5 

minimum  value  of  q£  -  6.4  x  10  W/m  deduced  for  the  center  of  the  circular 
region  to  reach  the  vaporization  temperature  is  larger  than  that 
calculated  earlier  because  the  assumed  higher  conductivity  of  the 
liquid  permits  faster  diffusion  of  heat  away  from  the  heat  source. 

Figure  4  shows  the  effect  of  different  thermal  conductiv¬ 
ity  values  on  the  location  of  the  liquid-solid  interface  when 
steady  state  prevails.  Note  that  a  higher  liquid  conductivity, 
while  the  conductivity  of  the  solid  remains  the  same,  results  in 
a  larger  metal  pool,  higher  £a  value,  if  the  center  of  the  circu¬ 
lar  region  is  to  reach  the  vaporization  temperature. 
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The  steady  state  temperature  distributions  of  a  stationary 
heat  source  are  of  interest  because  this  is  the  problem  geometry 
as  the  velocity  of  the  moving  heat  source  approaches  zero.  Fig¬ 
ure  5  shows  the  shape  and  location  of  several  isotherms,  includ¬ 
ing  the  liquid-solid  interface,  in  an  aluminum  substrate  for  two 
different  values  of  the  product  £a^.  These  are  steady 
state  isotherms  and  the  center  of  the  circular  region  has 
reached  maximum  temperatures  Tv  and  2130K  for  ^a  =  6.4  x  105 

5 

W/m  and  ga.  =  5  x  1 0  W/m,  respectively.  The  temperature  distri¬ 
butions  remain  the  same  in  these  dimensionless  plots  for  all 
values  of  £  and  a  as  long  as  the  product  c[a  is  kept  constant. 

Figure  6  shows  the  effect  of  increasing  the  product  £a  on 
the  steady  state  location  of  the  liquid-solid  interface.  The 
associated  maximum  temperature  at  the  center  of  the  circular 
region  is  listed  on  each  curve.  The  melt  pool  becomes  deeper 
and  hotter  as  the  product  of  the  uniform  absorbed  heat  flux  and 
the  radius  of  the  circular  region  increases.  Furthermore, 
ratio  of  the  melt  width  to  the  melt  depth  is  larger  than  one 
and  increases  with  decreasing  value  of  the  product  £_a. 

Figures  7  and  8  show  the  actual  temperature  distribu¬ 
tions  in  the  melt  pool  and  the  solid  along  the  z-axis  and  on  the 
z  =  0  plane,  respectively.  As  anticipated,  temperature  gradients 
along  the  z-axis  increase  with  increasing  values  of  the  product 
qa  and  decrease  with  increasing  distance  from  the  surface  of  sub¬ 
strate.  On  the  other  hand,  symmetry  on  the  z=0  plane  requires  a 
zero  temperature  gradient  with  respect  to  the  x  or  the  y  axis  at 
the  center  of  the  circular  region.  Note  that  the  temperature 
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gradients  at  the  edge  of  the  pool  on  this  plane  increase  with 
decreasing  values  of  the  product  This  information  will  be 

of  interest  when  relationships  between  cooling  rates  and  process 
variables  are  discussed  in  the  moving  heat  flux  problem  geometry. 

2 .  Quasi-Steady  State  Heat  Flow  — 

Moving  Heat  Flux 

A.  Uniform  Heat  Flux 

An  important  initial  finding  of  this  investigation  was  that 
the  three  variables,  absorbed  heat  flux  £,  the  radius  of  the  cir¬ 
cular  region  a^  and  the  velocity  of  the  moving  heat  flux  U,  could 
be  combined  into  two  independent  variables.  That  is,  the  dimen¬ 
sionless  temperature  distribution  in  the  liquid  metal  pool  and  the 
solid  substrate  remain  the  same  as  long  as  the  products  c[a  and  U£ 
or  U/q  are  kept  constant  while  the  individual  values  of  the  three 
variables  are  varied.  Consequently,  the  data  is  presented  herein 
in  a  general  form,  covering  a  large  range  of  process  parameters, 
in  terms  of  the  product  £a  or  U/q  and  the  dimensionless  variable 
Ua/2g. 

Figure  9  shows  the  dimensionless  temperature  distribution 
along  the  y-axis  for  different  values  of  Ua/2a.  For  small  val¬ 
ues  of  Ua/2a  ^0.003  the  heat  flux  is  moving  very  slowly  across 
the  substrate  in  the  y-direction  and  the  temperature  distribu¬ 
tions  are  almost,  but  not  exactly,  identical  to  those  shown  in 
Figures  5  to  8.  That  is,  there  is  little  distortion  of  the 
melt  pool  and  it  remains  almost  symmetrical  as  it  travels  across 
the  substrate.  On  the  other  hand,  increasing  the  dimensionless 
velocity  results  in  increasing  distortion  of  the  metal  pool  -  the 
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maximum  temperature  along  the  y-axis  shifts  toward  the  tail  end 
of  the  pool.  The  data  in  Figure  9  permits  determination  of  temp¬ 
erature  distribution  along  the  y-axis  for  a  wide  range  of  process 
variables.  For  example,  for  a  radius  of  the  circular  region 

а.  =  400pm  and  £a  =  6.4  x  108  W/m  the  maximum  temperature  reached 
at  steady  state,  JJ  =  0,  was  the  vaporization  temperature  of  alum¬ 
inum.  The  maximum  dimensionless  velocity  Ua/ 2a  ^1.0  in  Figure  9 
translates  into  an  actual  velocity  of  ^  0.42  m/s,  a  temperature 
T(0,0,0)  '\-2400K  and  a  maximum  temperature,  displaced  from  the  center 
of  the  circular  region,  Tmax  -  2540K.  Changing  the  product  cja_  to 

5 

4.1  x  10  W/m  while  Ua/2a  is  kept  constant,  results  in  a  signifi¬ 
cant  reduction  in  the  temperature  at  the  origin  T(0,0,0)  ^1 5 1  OK  — 
compare  this  with  the  steady  state  value  of  1730K  in  Figure  6. 

Examples  of  the  shape  and  location  of  several  isotherms, 
including  the  liquid-solid  interface,  for  given  values  of  the 
products  £a  and  Ua/2a  are  shown  in  Figures  10  and  11.  Figure  10 
shows  a  side  view,  plane  x  =  0,  while  Figure  11  shows  a  top  view, 
plane  z  =  0.  The  dimensionless  velocity  Ua/2g  =  0.75  would,  for 
example,  translate  into  actual  velocities  of  0.1  m/s  and  1  m/s  for 
radii  of  the  circular  region  of  ,vl  260ym  and  ,vl26pm,  respectively. 

The  corresponding  absorbed  uniform  heat  fluxes  that  result  in  £a  = 

б. 4  x  105  W/m  are  o  'v  5.1  x  108  W/m2  and  3  -v  5.1  x  109  W/m2, 
respectively.  Figures  10  and  11  also  show  significant  shifts  in 
the  geometry  of  the  isotherms  to  the  trailing  end  of  the  moving 
heat  source.  Also,  the  isothermal  surfaces  become  increasingly 
distorted  with  increasing  temperature  in  the  liquid  metal  pool.  A 
better  indication  of  these  observations  can  be  noted  in  Figure  12. 


This  Figure  shows  composite  views  of  the  isotherms  in  Figures  10 
and  n  along  with  side  and  top  views  of  the  isotherms  when  the 
product  ga^  is  the  same  but  the  heat  source  is  not  moving  -  the 
steady  state  temperature  distributions  when  no  further  melting  or 
solidification  occurs.  It  is  interesting  to  note  that  due  to  the 
high  conductivity  of  the  aluminum  substrate  the  distortions  in  the 
isotherms  are  not  nearly  as  pronounced  as  those  expected  in  lower 
conductivity  materials  such  as  iron  or  nickel.  This  point  was 
clearly  demonstrated  in  the  moving  point  source  cal cul ations  of 
Rosenthal  (3). 

Figures  13  to  16  show  the  effects  of  changing  the  variables 
qa  and  Ua/2g  on  the  geometry  and  location  of  the  liquid-solid  inter¬ 
face  -  the  liquid  pool.  First,  it  is  evident  that  there  is  little 
distortion  of  the  pool  at  low  velocities.  Second,  increasing 
the  product  Ua/2ot  results  in  a  corresponding  decrease  in  maximum 
pool  depth.  Third,  the  geometry  of  the  pool  is  more  speherical 
at  higher  products  of  £a  -  its  width  to  depth  ratio  increases 
with  decreasing  values  of  ga^.  Finally,  shallower  pool  geometries 
are  less  affected  by  changes  in  the  dimensionless  velocity. 

The  effect  of  changes  in  the  values  o.f  £a  and  Ua/2a  on  the 
cooling  rate  in  the  liquid  at  the  solid-liquid  interface,  G^'R*, 
and  the  variation  of  this  cooling  rate  along  the  different  axes  of 
the  cartesian  coordinates  are  shown  in  Figures  17  and  10. 

These  are  calculated  cooling  rates  during  solidification  of 
the  trailing  half  of  the  metal  pool.  Heating  and  cooling  rates 
from  the  point  of  view  of  a  stationary  observer  located  anywhere 

*6^  and  R  are  the  temperature  gradient  in  the  liquid  and  the  solid- 
liquid  interface  velocity  perpendicular  to  the  metal  pool  surface, 
respectively. 


in  the  heat  affected  zone  of  the  substrate  can  similarly  be  deter¬ 
mined.  The  cooling  rate,  G^ • R ,  can  be  alternatively  described  as 
the  product  of  U .  3  T  /  3  y ,  where  9  T  /  3  ,y  is  the  y-component  of  the  temp¬ 
erature  gradient  in  the  liquid  at  the  liquid-solid  interface.  The 
cooling  rate  is  a  maximum  along  the  y-axis  and  its  value  for  a 
given  radius  of  the  circular  region  increases  with  increasing  val¬ 
ues  of  the  ratio  U/q .  This  is  clearly  evident  in  the  plots  of  Fig- 

2 

ures  17  and  18  -  note  the  values  of  G|_-R  x  a  at  z/a  =  x/a  =  0. 

This  fact  can  also  be  deduced  from  the  steady-state  temperature 
distributions  in  Figure  8.  The  temperature  gradients  at  the  solid- 
liquid  interface  increase  with  decreasing  values  of  uniform  absorbed 
heat  flux  cj.. 

For  given  values  of  £a^  and  Ua/2ot  the  cooling  rate  continuously 
decreases  toward  the  edges  of  the  metal  pool.  Figures  17  and  18. 

This  is  expected  since  both  the  temperature  gradient  in  the  y-direc- 
tion  and  the  solid-liquid  interface  velocity  perpendicular  to  itself  reach 
a  finite  but  minimum  value  at  these  locations.  Examples  of  cooling 
rates  that  are  readily  calculated  from  these  Figures  are  as  follows. 
Assume  a  uniform  heat  flux  of  £  =  1.2  x  109  W/m2  absorbed  over  cir¬ 
cular  region  of  radius  a.  =  250ym  is  moving  with  a  velocity  IJ  0.5 
m/ s  >  M  ^  3  x  105  W/m  and  Ua/2a  -u  0.75.  The  calculated  cooling 
rates  from  Figures  17  and  18  at  y/a  =  0,  x/a  =  0.8  and  z/a  =  0.4 
are  9  x  10^K/sec,  n,  3.7  x  10^K/sec  and  n.  8  x  10^K/sec,  respec¬ 
tively. 

The  ratio  of  the  temperature  gradient  in  the  liquid  at  and 
perpendicular  to  the  solid-liquid  interface  divided  by  the  solid- 
liquid  interface  velocity  perpendicular  to  the  melt  pool  during 
»ol idi f ication ,  G^/R,  is  a  measure  of  the  stability  of  a  planar 
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interface  and  its  progressive  breakdown  into  cellular  and  dendri¬ 
tic  solidification  modes.  This  parameter  is  plotted  versus  the 
2/a  and  x/a  axis  for  different  values  of  the  variables  ga^  and 
Ua/2a  in  Figures  19  and  20,  respectively.  The  data  clearly  indi¬ 
cate  that  the  minimum  G^/R  value  consistently  occurs  along  the 
centerline  of  the  moving  pool  on  the  y/a  axis.  The  interface  vel¬ 
ocity  at  this  location  is  in  the  y-direction  and  assumes  its  maxi¬ 
mum  value  of  £.  The  temperature  gradient  in  this  location 
increases  with  decreasing  values  of  the  ratio  U/q ,  increasing  uni¬ 
form  absorbed  heat  flux,  for  reasons  already  discussed.  Therefore, 
the  ratio  G^/R  along  the  y-axis  increases  with  decreasing  ^  and 
increasing  velocity  U_. 

It  is  important  to  note  that  as  one  moves  along  the  back  of 
the  pool  on  the  y  =  0  plane  from  z/a  =  0  toward  the  bottom  of  the 
pool  the  interface  velocity  vector  both  rotates  and  changes  in 
magnitude  —  it  starts  out  at  its  maximum  value  of  U_  pointing  in 
the  positive  y-direction  and  continuously  decreases  to  its 
minimum,  small  but  finite,  value  at  the  maximum  pool  depth 
pointing  almost  in  the  negative  z-direction.  It  is  thus  clear 
why  the  G^/R  increases  with  increasing  distance  down  the  back 
of  the  pool,  increasing  z/a,  and  assumes  its  maximum  value  at  the 
maximum  pool  depth.  Figure  19. 

In  a  similar  manner,  simultaneous  rotation  and  decrease  in 
magnitude  of  the  interface  velocity  vector  occurs  in  the  z  =  0 
plane  as  one  moves  from  the  back  to  the  side  of  the  pool.  This 
explains  the  increasing  Gj_ / R  values  with  increasing  distance  along 
x/a  in  Figure  20. 
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Finally,  in  both  Figures  19  and  20  the  value  of  G^/R  decreases 
with  increasing  ratio  of  U/ q  at  constant  a^  or  increasing  £a  at  con¬ 
stant  lJa/2ot.  These  trends  are  readily  explained  by  following  the 
same  reasoni ng  as  that  presented  above  for  the  cooling  rate.  ! 

B.  Gaussian  Heat  Flux 

The  effect  of  changing  the  heat  flux  distribution  from  a  top 

hat  (uniform)  to  a  Gaussian  was  investigated.  As  previously  noted,  (2), 

if  the  total  absorbed  power  in  the  circular  region,  Q,  is  identical 

for  the  uniform  and  the  Gaussian  heat  flux  distributions,  then  the 

following  relationship  is  readily  deduced: 

q 

^uniform  ~  2.313  (48) 

where  qQ  is  the  absorbed  heat  flux  at  the  center  of  the  circular 
region  in  the  Gaussian  distribution. 

Examples  of  calculated  liquid-solid  interface  locations  for 
the  Gaussian  heat  flux  distributions  are  shown  in  Figures  21  and  22. 

5 

The  product  qQa/2 . 31 3  n,  4.65  x  10  W/m  resulted  in  a  maximum  steady 
state  T(0,0,0)  =  Ty  for  the  case  of  a  stationary  heat  source.  This 
value  is  lower  than  that  for  the  uniform  heat  flux  due  to  the  high 
concentration  of  absorbed  power  at  the  center  of  the  Gaussian  dis¬ 
tribution.  The  liquid-solid  interface  for  small  values  of  the 
dimensionless  velocity  Ua/2a  is  symmetrical  in  both  the  y  =  0  and 
z  =  0  planes.  Figures  21  and  22,  respectively.  For  a  given  temp¬ 
erature  in  the  center  of  the  circular  region  the  metal  pool  is 
shallower  for  the  Gaussian  heat  flux  distribution  -  compare  the 
solid  curve  for  Ua/2a  n,  0.03  in  Figure  21  with  the  curve  for 


c 

qa  =  6.4  x  10  for  the  same  dimensionless  velocity  in  Figure  13. 
This  is  expected  since  the  total  power  absorbed  in  the  Gaussian 
distribution  is  1 ower . 

Decreasing  the  product  qoa/2.313  results  in  colder  and  shal¬ 
lower  melt  pools  while  increasing  the  dimensionless  velocity 
shifts  the  trailing  end  of  the  pool  toward  the  negative  y-axis. 
These  observations  are  in  line  with  previous  findings  for  the 
case  of  uniform  heat  flux.  It  was  also  found  that  motion  of  the 
heat  source  has  a  more  pronounced  influence  on  the  melt  tempera¬ 
tures  than  in  the  previous  case.  For  example,  increasing  the 
dimensionless  interface  velocity  to  Ua/2a  'v  0.75  reduced  the  temp¬ 
erature  at  T(0,0,0)  from  Tv  to  ^1600K  for  the  value  of  q  a/2.313 
-x-  4.65  x  105  W/m. 

Finally,  general  trends  relating  cooling  rates  and  GL/R  values 
are  similar  to  those  previously  discussed  for  the  case  of  uniform 
absorbed  heat  flux  distribution. 

VI.  SUMMARY 

The  three-dimensional  temperature  distributions  in  the  melt 
pool  and  the  adjacent  heat  affected  zone  of  a  s^mi -inf  ini te  sub¬ 
strate  subjected  to  a  moving  directed  high  energy  source  can  be 
readily  determined  with  the  generalized  formulation  of  the  heat 
flow  equation  in  orthogonal  curvilinear  coordinates  coupled  to 
an  enthalpy  model.  While  numerical  computations  are  presented  for 
an  aluminum  substrate  subjected  to  moving  uniform  and  Gaussian 
heat  flux  distributions,  the  equations  developed  could  be  equally 
applicable  to  a  range  of  metallurgical  processes  previously  treated 
with  the  moving  point  source  equation.  It  is  shown  that  if  the  two 
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independent  variables  £a  and  Ua/2a  or  U/q  are  specified,  the  dimen¬ 
sionless  temperature  distributions  in  a  given  substrate  material 
remain  the  same.  Shortcomings  of  the  model  include  the  use  of  con¬ 
stant  thermophysical  properties  and  the  fact  that  convection  in  the 
metal  pool  is  not  taken  into  consideration  except  by  arbitrarily 
increasing  its  thermal  conductivity.  On  the  other  hand,  the  gen¬ 
eral  trends  that  can  be  deduced  for  a  given  energy  source  and  sub¬ 
strate  material  should  permit  a  more  systematic  approach  to  the 
variation  and  control  of  the  process  variables  in  order  to  achieve 
the  desired  heat  flow  conditions  during  melting  and  solidification. 
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APPENDIX 


An  alternate  method  to  the  general  formulation  of  the  finite 
difference  representation  of  the  heat  conduction  equation  in  moving 
orthogonal  curvilinear  coordinate  system,  equation  (1),  is  con¬ 
sidered  here.  In  this  method  an  energy  balance  is  applied  to  a 
control  volume  of  finite  size  in  oblate  spheroidal  coordinates 
moving  with  velocities  v^,  v^  and  v^,  see  Figure  2.  The  resultant 
energy  balance  is  approximated  using  finite  differences  and  is 
shown  to  converge  to  an  identical  formulation  as  that  given  in 
equations  (19)  and  (20).  This  approach  is  used  to  both  verify  the 
general  formulation  of  equation  (1)  developed  in  this  investigation 
and  to  render  improved  physical  interpretation  of  the  various  coef¬ 
ficients  given  in  equation  (20). 

The  moving  control  volume  element  about  a  point  (i,j,k)  is 
illustrated  in  Figure  2.  The  energy  balance  for  this  volume  element 
having  a  total  source  strength  of  PAV  is  carried  out  by  considering 
the  total  rate  of  heat  transfer  through  each  surface  by  conduction 
and  due  to  the  motion  of  the  volume  element. 

Let  Q  and  Qv  denote  the  rate  of  heat  transfer  entering  or 
leaving  the  volume  element  by  conduction  and  by  motion  of  the  vol¬ 
ume  element,  respectively.  The  energy  balance  for  the  quasi -steady 
state  case  under  consideration  is  given  by: 


(Qi-l/2,j,k  “  Qi  +  l/2,j,k^  +  ^i,j-l/2,k 
+  ( Q i , j , k - 1 / 2  “  Qi,j,k+l/2)  +  (qi+l/2,j,k 
+  {Cri,j  +  l/2,k  "  Qi,j-l/2,k)  +  (Qi,j,k+l/2 


Qi,j+l/2,k) 

Qi-l/2,j,k) 

q'i,j,k-l/2} 


+  PAV  =  0 


(Al) 


t 


In  what  follows  the  values  of  Q  and  Qv  for  the  1-1/2  and  1+1/2 


faces  are  evaluated  as  examples.  All  the  other  terms  in  equation 
(Al)  can  then  be  similarly  obtained.  It  is  shown  that  the  sum  of 
all  these  terms,  after  appropriate  manipulations,  results  in  equa¬ 
tion  (19)  with  the  coefficients  defined  in  expression  (20). 


Qi-l/2,j,k  =  “k  -Ti,J,k-" 

(A  Vi-1/2 


AAi - 1 /2 


(A2) 


Appropriate  substitutions  of  equations  (5)  to  (7)  and  (16)  into 
equation  (A2)  and  some  manipulation  gives: 


Qi-l/2,j,k  =  '  Cl(Ti,j,k  "  Ti-l,jV 

Similarly, 

Qi  +  l/2,j,k  =  ~  Wl.j.k  “  Ti,jV 


akcoshn  .sin?;  .AnA£A<j>  (A3) 

■  J 

akcoshq.:  sin£.AnA?;A<i>  (A4) 

*  J 


The  coefficient  and  C2  are  those  defined  in  expressions  (20). 
Similar  expressions  to  (A3)  and  (A4)  are  readily  developed  for  the 
rates  of  heat  transfer  through  the  other  four  faces  of  the  volume 
element  due  to  conduction. 

The  rate  of  heat  transfer  entering  the  i-1/2  face  due  to  the 
motion  of  the  volume  element  is: 


AA, 


Qi-l/2,j,k  =  pCpTi-l/2,j  ,k  Vr!._1/2  “"i-1/2 


(A5) 


Substitution  of  equations  (6),  (16)  and  (18)  into  equation 
( A  5 )  gives: 


Qi 


i-1/2, j,k 


aUpC  s  i  n£  .  s i no . 

_ _  JtL _ _ E 

kAn 


C  T  i  - 1  /  2  ,  j  ,  k 


,  J  • 


sinhni_1/2coshni_i/2 


coshn^ 


•  akcoshn^sin^ .AnA^Ap 


(A6) 


I 


V 
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aUpC  s i n£ i s i n<j>. 

«)-l/2,j,k  ■  -  p  kA'n -  <Ti-l/2,j.k>  •  •kc«shni«ln£j4n»?A» 


cosh2n_-  . 

[5in,"'i  -  -coshn,  “r] 


Similarly: 


1  i  + 1  /  2 ,  j  ,  k 


aUpCpsing  ^  s  i  n  <f>  k 


( T i  + 1  / 2 ,  j  ,  k )  *  akcoshn^in?  .AnACA<!> 


cosh2n,-  * 

[sinh,1i  +  ~coshnT~  -T] 


Subtraction  of  equation  (A7)  from  (A8)  gives: 


Qi+l,j,k  "  Qi-l/2,j,k 


aUpCpsingjSin({)k  »  akcoshn^ si ngjAnA£A<{> 


sinhri.  cosh2n 

2Arj  ^Ti+l,j,k  "  Ti-l,j,k^  +  coshn* 


r(T<^] 


Substituting  eauat-ions  (A3),  (A4),  (A9)  and  similar  terms  for 
the  other  four  faces  into  equation  (Al)  and  dividing  both  sides  by 

akcoshn-s  in?  .AnA£;A$ 

•  J 

gives  an  identical  equation  to  expression  (19).  Note  that  the  sum 
of  the  second  terms  inside  the  brackets  of  equation  (A9),  the  coef¬ 
ficient  to  T,  .  becomes  zero  in  equation  (Al).  This  essentially 
implies  that  the  sum  of  the  combination  of  area  and  velocity  terms 
multiplied  and  differentiated  with  respect  to  each  axis  is  zero  or 


pC  hv 


n  =  1  3ul 


-L)  «  0 


(A10) 


which  was  the  a  priori  assumption  made  in  the  derivation  of  equa¬ 
tion  ( 1 ) . 
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a 


A 

C 


x  ,y  ,z 


NOMENCLATURE 

radius  of  the  circular  region,  m  or  ym 
area  of  the  element 
integration  constant 
specific  heat,  JKg'^K”^ 
fraction  liquid 

temperature  gradient  in  the  liquid  at  the  liquid- 

solid  interface,  Km”^ 

scalar  factor 

specific  enthalpy,  JKg-1 

heat  of  fusion,  JKg-^ 

thermal  conductivity,  Jm~^s"^K~^ 

_  3 

rate  of  heat  generation  per  unit  volume,  Wm 

_  o 

absorbed  heat  flux,  Wm 
rate  of  total  absorbed  heat,  W 
interface  velocity,  m  sec-1 
arc  length 
time,  s 

temperature,  K 

ambient  temperature,  K 

melting  temperature,  K 

vaporization  temperature,  K 

coordinate  axis 

velocity  of  heat  source,  ms~^ 

volume,  m^ 

velocity,  ms”^ 

cartesian  coordinates 


2  - 1 

a  thermal  diffusivity  ( k/pCp) ,  m  s 

n,C»4>  oblate  spheroidal  coordinates 

6  dimensionless  temperature  variable 

ip  dimensionless  enthalpy  variable 

_  3 

p  density.  Kg  m 

Subscripts 

i,j,k  nodal  point  subscripts  in  n  and  £  and  4»  directions, 
respectively 
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TABLE  I 


PROPERTIES  OF  THE 
ALUMINUM  SUBSTRATE 


=  1067 


0  Kg”^  K”\  Specific  Heat 


=  3.95  x  10' 


=  2723 


=  2545 


0  Kg“  ,  Latent  Heat  of  Fusion 

Wm"1  K'\  Thermal  Conductivity 
K,  Melting  Temperature 

K,  Vaporizati  n  Temperature 
-  3 

Kg  m  Density 


c  o  _  i 

8.4  x  10  m  sec  ,  Thermal  Diffusivity 


*  Averaged  from  298K  to  the  Vaporization  Temperature 
**  Averaged  from  298K  to  Melting  Temperature 


Figure  1 . 


Schematic  illustration  of  laser  beam  -  substrate 
geometry  during  rapid  surface  melting  and  solid¬ 
ification. 
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Figure  2 . 


Yhree-d ircens iona 1  moving 
dinate  system  and  problem 
critized  space  domain. 


oblate  spheroidal  coor- 
qeometry  showing  a  dis- 


10"2  10”1  10°  101  10 


10° 


10‘J 


a/2/a$t 


Figure  3.  Temperature  at  the  center  of  the  liquid  zone  of 
a  semi-infinite  solid  substrate  during  surface 
melting  as  a  function  of  uniform  absorbed  heat 
flux,  radius  of  the  circular  region  and  time. 


k[T(0,0,0) 


J 


FRACTIONAL  DISTANCE,  z/a 


Figure  5.  Steady  state  location  of  several  isotherms, 

including  the  liquid-solid  interface  (T=933K) 
for  two  different  products  of  qa . 


FRACTIONAL  DISTANCE,  z/a 


X 


FRACTIONAL  DISTANCE,  y/a  or  x/a 


-2.0  -1.6  -1.2  -0.8  -0.4  0  0.4 


0.8  1.2  1.6  2.0 


Figure  6.  The  effect  of  the  variable  £a^  on  the  steady 
state  geometry  and  location  of  the  liquid- 
solid  interface.  The  maximum  temperature  at 
the  center  of  the  circular  region  associated 
with  each  curve  is  also  listed. 


TEMPERATURE, 


TEMPERATURE 
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FRACTIONAL  DISTANCE,  y/a  or  x/a 


Figure  8.  Steady  state  radial  temperature  distributions 

in  the  melt  pool  on  the  z=0  plane  for  different 
values  of  the  variable  qa . 


to°Z/n*(  3/  HV 


-2.4  -1.6  -0.8  0  0.8  1.6  2.4 

FRACTIONAL  DISTANCE,  y/a 


Figure  9.  Quasi -steady  state  dimensionless  temperature 
distributions  along  the  y-axis  of  a  moving 
uniform  heat  flux  £  absorbed  over  a  circular 
region  of  radius  a.  The  velocity  of  the  heat 
source,  JJ_,  is  in  the  positive  y-direction. 


A  side  view,  x  =  0  plane,  of  the  location 
and  geometry  of  several  isotherms,  includ 
ing  the  liquid-solid  interface. 
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Figure  13.  A  side  view,  x  =  0  plane,  showing  the  effects 
of  changing  the  values  of  £a  and  Ua/2a  on  the 
shape  and  size  of  the  molten  region. 


FRACTIONAL  DISTANCE,  x/a 


FRACTIONAL  DISTANCE,  y/a 
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Figure  14.  A  top  view,  z  =  0  plane,  showing  the  effects 
of  changing  the  values  of  £a  and  Ua/ 2a  on  the 
shape  and  size  of  the  molten  region. 


FRACTIONAL  DISTANCE,  z/a 


Figure  15.  A  side  view,  x  =  0  plane,  showing  the  effects  of 
changing  the  values  of  ^a  and  Ua / 2a  on  the  shape 
and  size  of  the  molten  region. 
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a  =  0  75 


UNIFORM  HEAT  FLUX 
qa  =  3  x  10s  W/m 
qa  =  6.4  x  105  W/rn 
QUASI -STEADY  STATE 
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FRACTIONAL 


0.8  1. 
DISTANCE,  z/a 


Variation  of  the  product  of  cooling  rate  at 
the  solid-liquid  interface  and  a_^  with  melt 
depth  along  the  trailing  half  of  the  pool  in 
the  y  =  0  plane  for  different  values  of  the 
independent  process  variables. 


53 


\ 


0  0.2  0.4  0.6  0.8  I. O'  1.2  1.4  1.6 
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Figure  78. 


Var  at;  on  of  the  product  of  cooling  rate  at  the 

?  '  qU1Vn1TfaCe  and  a"  with  ne1t  width 

a  ong  the  side  of  the  pool  in  the  z  =  0  plane  for 

abl  es!'ent  Va  ues  of  the  ^dependent  process  vari- 


Figure  19.  Variation  of  the  ratio  Gi/R  at  the  solid-liquid 

interface  with  melt  depth  along  the  trailing  half 
of  the  pool  in  the  y  =  0  plane  for  different  val¬ 
ues  of  the  independent  process  variables. 


FRACTIONAL  DISTANCE,  x/a 


Figure  20.  Variation  of  the  ratio  GL/R  at  the  solid-liquid 
interface  with  the  melt  width  along  the  side  of 
the  pool  in  the  z  =  0  plane  for  different  values 
of  the  independent  process  variables. 
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Figure  21.  A  side  view,  x  =  0  plane,  showing  the  effects  of 
changing  the  values  of  £a  and  Ua/ 2a  on  the  shape 
and  size  of  the  molten  region  of  an  aluminum  sub¬ 
strate  subjected  to  a  Gaussian  heat  flux  moving 
with  constant  velocity  U  in  the  positive  y-direc- 
tion. 
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